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INTRODUCTION 


JHU/APL acts as technical direction agent 
for US Navy weapon system programs 


A key task of APL’s Guidance, Control, and 
Navigation Systems Group is the evaluation 
or conceptual design of missile guidance and 
control systems 


Analysis and design work requires a flexible. 
interactive linear modeling tool 


PC-MATLAB resident on 80386 engineering work- 
stations provides such a tool 


Work presented here shows general attributes of 
MATLAB, demonstrating use of PC-MATLAB/386 
for linear design of a flexible missile autopilot 
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MATLAB BACKGROUND 


MATLAB (MATrix LABoratory) provides an interactive, 
matrix-oriented environment 


■ MATLAB is based on the EISPACK and UNPACK routines 
for matrix computations 


■ PC-MATLAB/386 is a high-performance MATLAB 
implementation for 80386-based workstations 


MATLAB built-in functions, plus higher-level 
functions developed for control system calcula- 
tions, allow for effective controls design 
studies 


144 


HARDWARE AND SOFTWARE CONFIGURATION 



■ COMPAQ 386/20 computer 

■ Weitek 1167 numeric coprocessor 

■ PC-MATLAB/386 with Control Systems Toolbox 
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Interactive, high-level command environment 


Very high processing speed 


Easy extensibility via user-defined functions 
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A MATLAB INTERACTIVE COMMAND LINE EXAMPLE 


» k = lqr(a,b,q,rho*r); eig(a-b*k), y = step(a-b*k,b,c,d,l,t); plot(t r y); 


■ The single line above, typed at the MATLAB 
command line prompt, does several things: 

- Computes a quadratic regulator gain vector 

- Displays the closed-loop eigenvalues — often 

useful for confirming that actuator band- 
width requirements are not excessive 

- Computes and plots a unit step response 


■ By varying the control cost (rho) above, a very 
large family of compensators may quickly be 
considered 


■ The above command line suggests the power and 
utility available from a high-level, inter- 
active matrix language 
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PC-MATLAB/386 PROCESSING SPEED 


■ MATLAB's UNPACK Benchmark: 460 double precision KFLOPS 


This processing speed is: 

- 25 x faster than standard PC/AT 

- 6 x faster than Mac II 

- 3 x faster than MicroVax II 


■ Implication: the fast response time resulting from 
such performance allows for truly interactive 
design iterations on complex control laws 
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MATLAB EXTENSIBILITY 


■ User-defined functions may be developed 
through creation of simple text files 


■ Some typical user-defined functions: 

- Frequency-response plotting routines 

- Application-specific linear transformations 

- Multivariable Nyquist criterion 

■ Complex state-space or transfer-function 

models also defined through user text files 
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AN EXAMPLE OF A USER-DEFINED COMMAND FILE 


■ Below command set calculates and plots the 
maximum and minimum singular values of 
a plant and observer-based compensator, 
for a loop broken at plant input 


function [smin.smax] = avdinput(a ,b, c, kcon,kobs, w) ; 

X 

jay = aqrt(-l) ; 

[nn, xx]=sise(a) ; i2=eye(nn); [n*,xx] =sire(c*a*b) ; phi = *(**i2-a)’; 
for i = 1 : nc ; 

s = w(i)*jay; phieval = eval(phi); 

ga = c/phleval*b ; ks = kcon / { phieval+b*kcon«-koba*c ) * koba; 
xx=avd(ks*ga ) ; amin( i )=xx(ng) ; amax( i )=xx( 1 ) ; 
end ; 

X 

X 

X convert to decibels and plot output 
X 

amin=20*logl0( amin) ; amax=20*logl0 ( anax ) ; 
semilogx( w, amin, w, anax, ' r— ' ) i grid ; 

title(‘Max and Min Singular Values; Loop Broken at Plant Input *); 
xlabel (* Frequency (rad/aec)'); ylabel ( 'Magnitude (db)*); 


Procedure requires only eleven lines of 
executable MATLAB code 




CLASSICAL CONTROL CAPABILITIES 


Frequency response 


Root locus 


Nyquist plots 


Development of dynamic compensators 
(lead-lag, notch filters, etc) 
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MODERN CONTROL DESIGN EXAMPLE 


Design plant describes tactical missile at 
a high-altitude flight condition 


Design plant includes single-plane rigid- 
body dynamics and effect of first flexible 
mode on sensed pitch rate 


Objective is to develop an autopilot to track 
commanded accelerations 


Design challenge is to achieve high closed-loop 
bandwidth in presence of low-frequency 
bending modes 
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Establish design goals for closed-loop 
responsiveness and stability 


Develop full-state feedback (LQR) gains 
for design plant 


Define linear observer to reconstruct full state 
vector 

- Use "robust observer" design (Doyle and Stein. 

1979 IEEE Transactions on Automatic Control) 

- Adjust observer gains to recover original LQR 

loop transfer in desired frequency range 
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DESIGN PLANT MODEL 



Fifth-order state vector x; x = Ax + feu 
- = far a / s <lf/s q f ] 


First three state variables are associated with 
rigid-body airframe; the last two describe 
flexible mode dynamics 


Rate gyro measurement: [1 0 0 0 1) • x 


(Integrated) accelerometer measurement: [0 0 1 0 0] * x 


1 . 0000e+00 
0 
0 
0 


~2 . 3557©+02 
0 

2 . 6l58e+00 
0 
0 


1 . 7967e+02 
0 

-1 . 995l©+00 
0 
0 


0 

0 

0 

0 

-2.4649e+04 


1 . 0000«+00 
“3. 1400e-f00 


“2 . 803le+02 
0 

9.2587e-02 
. 0 
3 . 0723e+02 


154 


O O O 


SOME OBSERVATIONS ON DESIGN PLANT MODEL 


■ Feedback of the first three states describes 
a very standard (rigid-body) autopilot 
topology, used by tactical missiles since 


1950’s 


■ Open-loop plant is characterized by lightly 

damped airframe (weathercock) poles, and 
by bending mode poles 

- Airframe pole frequency lies at 

nominal 2.5 Hz 

- Bending mode has nominal 25 Hz natural 

frequency 


Desired autopilot crossover frequency here 
will lie near the bending mode frequency 


Pitch Rate (deg/sec) 


EFFECT OF STRUCTURAL MODE ON SENSED PITCH RATE 
(RATE GYRO MEASUREMENT) 




CONTROLLABILITY AND OBSERVABILITY PROPERTIES OF PLANT 


■ System (A,b) is controllable 


■ System is unobservable if rate gyro alone, 
or accelerometer alone, is used as the 
measurement to reconstruct state vector 


■ Both sensor outputs thus should be used in the 
observer design 


Approach taken for this application: 

- Define a (non-square) design plant having 

one input (fin deflection) and two inde- 
pendent outputs (gyro and accelerometer) 

- Use extensions of loop transfer recovery 

(Williams and Madiwale, 1985 ACC) valid 
for non-square systems 
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FREQUENCY RESPONSE OF FULL-STATE FEEDBACK (LQR) SYSTEM 
(LOOP BROKEN AT PLANT INPUT) 



OBSERVATIONS ON LOOP TRANSFER RECOVERY PROCEDURE 


■ For this application, recovery at both the 
(rigid-body) airframe and bending mode 
frequencies may only be achieved with very 
high observer gains 


■ For practical ranges of observer gains, recovery 
at airframe frequencies is obtained at the cost 
of lessened robustness in the structural mode 
frequency range 


■ Use of a set of user-defined MATLAB files, to 
implement a range of observer gain calcu- 
lations, makes evaluation of this robustness 
tradeoff straightforward 


IS 9 


RECOVERY OF DESIRED FULL-STATE FEEDBACK SYSTEM 
WITH MODEL-BASED COMPENSATOR 



Achieved Acceleration (Gees) 


ACCELERATION STEP RESPONSE OF FINAL COMPENSATOR DESIGN 
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Deg/sec 


RESPONSE OF FLEXIBLE MODE STATE DURING 
ACCELERATION STEP RESPONSE 


Pitch Rate Response Due to Flexible Mode 
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Arrei PRATION STEP RESPONSE FOR CASE WHEN 
BEND^G MODeTs PERTURBED TO 25 % LOWER VALUE 







Deg/sec 




SUMMARY OF DESIGN RESULTS 


■ Model-based compensator yields a high-bandwidth 
autopilot, which is robust to at least a 25% 
perturbation in bending mode frequency 


■ A number of issues still not addressed: 

- Detailed noise sensitivity assessment 

- Effect of higher-frequency structural modes 

- Phase lag from actuator dynamics 

- Effect of structural modes on accelerometer 

measurement 

- Tolerance to uncertainties in aerodynamics 


Above concerns could also be addressed using MATLAB 


SUMMARY: MATLAB APPLICABILITY FOR 
CONTROL DESIGN OF FLEXIBLE SYSTEMS 


MATLAB provides the necessary tools for a 
variety of control system design techniques 


Extensibility of MATLAB allows development 
of tools to implement recent modern control 
design methods, including loop transfer 
recovery 


■ Implementation for 80386-based machines (PC- 
MATLAB/386) has very high performance, 
allowing for interactive control design of 
complex systems such as flexible structures 


Any flexible structures control problem which 
can be cast into a state-space framework 
may benefit from design work with MATLAB 
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